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1, INTRODUCTION 

As one of the alternative of the renewable energy technology, wind turbine is now not only the most 
economic, but also benefit us socially and environmentally. The global generation of wind energy increased 
up to 435 GW in 2017, which is nearly 10% of overall world power generation [1]. In EU, wind energy has 
produced more capacity than coal power and accounts for nearly one fifth of installed power generation 
capacity in 2016 [2]. The result of global community pushing towards better environment through better 
ways of producing energy has produced many technological breakthrough that not only increase more 
renewable energy production, but also lower its generation cost. 

However, there are still challenge that need to be met to fully utilised wind energy potential. Chief 
amongst them is the fact that wind energy varies according to the speed of the wind. The variability of the 
grid connected wind power generation can be absorbed by the conventional power generators, but as the 
penetration of the wind capacity increased, rigorous studies and research must be done to anticipate and find 
solutions to the problem that might occurs [3]. 

The fluctuation of the wind generation could increase the temperature of the existing thermal 
generators beyond normal operating temperature [4]. As the wind generation varies throughout the time, the 
conventional generators need to act to compensate to meet the energy demand. This action will increase the 
temperature of these machines. The effect increased with the increased penetration of wind energy in 
the grid. 

In this paper, simulation study in a microgrid will be done covering 24 hours period which is then 
repeated on different set of wind speed data. Four different sets of wind speed according to historical data 
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from four different times of the year are used. The results from the simulation are then analysed and 
discussed in order to investigate the thermal effect of the wind energy fluctuations to the conventional 
generator. A microgrid where renewable energy is connected are more vulnerable to the intermittent nature of 
wind generation as the burden of compensating it will be largely fall upon the conventional generators in 
the system. 


2. WIND ENERGY 
2.1. Wind Speed 

One major characteristics of wind in any parts of the world is the changing patterns of its speed in 
different time of the year, or month, or season or day [5]. The speed is also varies according to different parts 
of the world geographically [6]. Although seasonal variations is more predictable than long term, year or 
decades variations, ther is also variations within a day or diurnal variations [7]. Shorter variations in minutes 
or even seconds, turbulences, also must be considered in studying the design and performance of a wind 
farm. Therefore in order to be able to predict or estimate wind energy production, it is important that this 
phenomenon be analysed and understood. 

As the wind speed changes according to the nature, the wind energy generation follows according to 
the speed change, instead of reacting to the energy demand as do the conventional generators. Figure 1 shows 
the 10 years average speed of wind in UK from 2002-2011[8] . In this example, the wind speed during the 
summer is noticeably slower than it is during the cold winter season. Although monthly pattern generally 
follows the trend of energy demand in UK where higher demands of energy are in the winter months, the 
wind speed patterns do not mirror the daily electrical demand [9]. 
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Figure |. Ten years average UK wind speed from 2002-2011 and 2016 


Although year to year variations are hard to predict and not well understood, the seasonal yearly 
variation can be presented in terms of probability distribution [10]. Wind speed characteristics during a 
period of time is shown to be following certain probability distribution called Weibull probability distribution 
[11]. It 1s govern by the following equation: 


k (v\k-1 =(-)* 
f(v,4,k) = AG ee ee (1) 
0,v <0 


where v is the wind speed, k is the shape parameter that gives form of the distribution; usually lies between 2 
and 3, and A 1s the scale parameter which indicates the width and peak location of the distribution. Rather 
than observing the changes of the wind speed by its mean value, the Weibull distribution let us to examine its 
whole characteristics [5]. As shown later in table 1, the values of k and A varies throughout the year, ranging 
from 2.08 to 2.53 for k and 11.52 knots to 16.25 for A. 


2.2. Wind Energy Generation 
The power output of a given wind turbine can be expressed by the following equation: 


1 

Po=5CpAy, B)patReV° (2) 
where P, is the power harvested from the wind, C, 1s the power coefficient of the turbine, Ay is the tip speed 
ratio, B is the blade pitch angle, R; is the wind turbine radius,V is the wind speed and fp, is the air density. The 
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power coefficient is a function of Ay and B [12]. C, is the proportion of amount of wind energy can that 
extracted, which has upper limit of 59% theoretically. Currently, modern wind turbines have values between 
45% to 50% [13]. As the power produced is directly proportional to the cube of the wind speed, any change 
in the speed will greatly affected the power output of any wind generator [14]. 


3. RESEARCH METHOD 
3.1. Thermal Effect Calculation 

In any generators, heating effect must always be taken into consideration as it will not only 
determine the operation of the generators but also the lifespan of the machine. 

The heating effect of a generator comes from the excitation current If and also stator current I. 
According to [15][ 16], the thermal energy W is given by the following equations: 


TU 
W = 1,B2 = 1,NfE? + 1,N212X4 — 2T,N,NEIXqcos (5 —§- 9) (3) 


where Tk is the proportionality constant relating W and Be2 , Be 1s the resultant flux density, Is 1s the stator 
current, Na and Nf are number of turns in the windings, Xd is the synchronous reactance, E is the no-load 
voltage, 5 is the rotor angle and @ 1s the power factor. 

The equation can be re-arranged to give: 


W = [q]'[A] (4) 
where [A]T =[al a2 a3], [q/T=[ql q2 q3] 
and ay = T,,N?, az = T,.N2,a3 = NaNel, 
TU 
q, = E2,q> = 12X2,q3 = —2EIXqcos (5 -5-@) 
As the heat produced by the generator increase its temperature, it is dissipated through the 
surrounding. If the rise of temperature between operating condition T and no load temperature Ts is T-Ts, 


and a is the transfer index for temperature and heating, therefore the following equation can be used to find 
the rate of increasing temperature. 


H = W(t)- a(T -T,) (5) 


H is the thermal inertia constant. The typical value of T; 1s 40°C. As the temperature (dy/dt) does not change 
in steady-state condition, then 


when t = 0,7) = me + Ts; when t = 0,7,, = me + T, 
Substituting To and T. into (5), then 
W(t) W(t) @). 
Se (“SY -1,+7,) eH + T, (6) 


The equation is used in the simulation to calculate the temperature of the machine where a = 1/300 
MW/°C, H/a = 600 s, H = 2 MW.s/°C, Xd ~ Xad = 1.305 pu, [A] = [9.1a, 33.20, 17.40] [15]. 


3.2 Cooling System Calculation 

All conventional thermal generators need to dissipate the heat generated by its machine core. This is 
the function of its cooling system. Main methods of cooling are once through systems, wet recirculating or 
closed-loop systems and dry cooling. 

The simplest method is once through systems where water from nearby sources is circulated in the 
condensers to absorb heat and then discharged locally. Wet recirculating systems recycle the water by using 
cooling towers to cool the warm water and therefore use less water. Dry cooling systems use air instead of 
water to cool the machine. There is also hydrogen cooling, where hydrogen is used to absorb the heat. In this 
paper, water cooling systems are modelled to be simulated as cooling system of the generators, which has 
better heat removal capability. 
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The heat absorbed by the cooling water flowing through the length dx of the generators is given by 
the following expression [4]: 


dQ, = a(S,(T — T,,)dx, dQ, = WC,dT,, (7) 
where S; is the perimeter of transversal section of waterway, W is the cooling water mass rate, Cy is the 


specific heat of cooling water at constant pressure and dTy is the difference of water temperature of length dx. 
Also: 


a, = 0.023Re°® prot (8) 
Dy 


where a 1s the convection coefficient of water, K,, is thermal conductivity of water; Dy is the equivalent hyd 
raulic diameter of waterway; Re is Reynolds number, Pr is Prandtl number. Both respectively are given by: 





HA Kw 
At steady state, 
@.51(T —Ty)dx =WCydTy, = “ax, as “88 = NTU (10) 


As both sides of the equation integrated over Ty and x, we get: 


NTU 


Ty =T—(T—Timje te (11) 
When x = Le, Tw = Tout, so: 
Laie sai (aos (T , Ewes: (12) 


In the simulation, data from [4] are used where W = 25.5kg/s, Cy = 4.1813x10°kg.K/J, Kw = 0.6 
W/m/K, p = 8.01 x 107 Pa-s, Tin = 45°C, Le = 12m. 


3.2. Simulation 

The simulation was done using Matlab and Simulink software. The microgrid consists of a 1SMW 
conventional synchronous generator and 4.5MW windfarm connected in parallel to a load through 
transmission line. The load is residential of IOMW, 0.95 pf. The generator acts as the base generator. The 
wind farm model is a simplified model that follows linear relationship with the wind. It will trip from the grid 
when the wind speed go beyond maximum permissible wind value, and reconnect when it go below it. The 
diagram of the system is shown in Figure 2. 





Generator Wind Farm 
15 MW 4.5 MW 





Load 10MW 
Pf =0.95 














Figure 2. Microgrid sytem in the simulation 


Wind speed data from four different days of the year corresponding to different seasons of the year 
in UK are used as shown in Figure 3 [17]. The value for the shape parameter, k and the scale parameter A are 
different for different seasons of the year as shown in Table 1 [18]. These values are used to generate wind 
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speed at 10 seconds interval according to Weibull distribution function. The simulations were run 
continuously for twenty four hours with constant load. The power generated by the windfarm and 


conventional synchronous generator were observed and the temperature of the machine and its water cooling 
were calculated. 
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Figure 3. 24 hours hourly wind speed data for Plymouth 


Table 1. Values of k AND A for different seasons at Plymouth Weather Station 


Season k X 
Winter (DJF) om | 16.25 
Spring (MAM) 2.4 13.83 
Summer (JJA) 258 11.52 
Autumn (SON) 2.08 13.73 


4. RESULTS AND ANALYSIS 


As the wind power generation varies according to the variability of the wind speed, the synchronous 
generator adjust its power output to meet the demand. The graph in Figure 4 shows both generators output 
during the first two hours of the simulation during the winter day. 
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Figure 4. Wind generation output and conventional generator output at 00:00 to 02:00 


The results of the machine temperature rise and its coolant are shown in Figure 5. The graph shows 
the fluctuation of the temperature of the machine in the period of 24 hours at different times of the year. As 
the wind speed and wind generation in the system fluctuate, the synchronous generator need to compensate 
the total generation to meet the demand. The machine temperature increased (and decreased) as a result of 
this. 

This results are then compared to the machine temperature when there is no wind generation in 
connection. Although the generator has to meet all the demand by itself, hence works closer to its full 
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capacity, there is no significant temperature rise in the core of the machine in comparison with all the other, 
which remains stable at around 50°C throughout the simulation. 
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Figure 5. Machine and cooling water temperature at different time of the year 


According to the calculation, the mean generator temperature for the day in winter, spring, summer 
and autumn is 80.5°C, 74.1°C, 60.7°C and 73.4°C. The higher temperature in the winter is due to more 
changes and fluctuation of wind speed during the winter day. It can be seen from the graph that the increase 
of the machine temperature can increase even though it 1s generating less power than usual. For example, the 
average temperature of the conventional generator during winter day is more than during the summer day, 
although more power can be generated by the wind farm during winter. 


5. CONCLUSION 

The paper discuss the effect of wind power fluctuation to the heating of conventional generator in a 
microgrid. The temperature increase in the machine and the cooling water can be analysed after simulation of 
the microgrid 1s carried out using the method discussed. 

The result of the simulation shows that there will be increase in the machine temperature as the wind 
generation fluctuate according to the nature of the wind speed. The rate of temperature increase also differs 
according to different values of k and A of the wind speed, which differs according to different time of the 
year in four season weather. 

Further studies need to be done in a bigger grid system to investigate this effect at different level of 
wind energy penetration. Also, a variable load corresponding to variable energy demand can also be 
incorporated in the future simulation to analyse the combine effect of the changes of the wind power 
generation and the load demand. 
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